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Abstract 



The Fourier transform is often used to connect the Lorentzian energy distri- 
bution for resonance scattering to the exponential time dependence for decaying 
states. However, to apply the Fourier transform, one has to bend the rules of 
standard quantum mechanics; the Lorentzian energy distribution must be ex- 
tended to the full real axis — oo < E < oo instead of being bounded from below 
< E < oo ("Fermi's approximation"). Then the Fourier transform of the ex- 
tended Lorentzian becomes the exponential, but only for times t > 0, a time 
asymmetry which is in conflict with the unitary group time evolution of stan- 
dard quantum mechanics. Extending the Fourier transform from distributions 
to generalized vectors, we are led to Gamow kets, which possess a Lorentzian 
energy distribution with — oo < E < oo and have exponential time evolution 
for t > t = only. This leads to probability predictions that do not violate 
causality. 
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PACS number(s): 03.65.Db, 03.65. Ca, 32.70. Jz, 06.30.Ft 



1 Introduction 

In this paper we would like to draw a connection between two theoretical questions concerning 
the decay of quasistable states. 

One question concerns the relation between the two experimentally-independent signa- 
tures of quasistable states, the Breit-Wigner or Lorentzian lineshape of resonances as mea- 
sured in the cross section and the exponential decay rate measured for decaying states. 
Explicitly, how can the standard relation T = h/r (based on the Weisskopf-Wigner approx- 
imation JTJ), which connects the lineshape parameter T to the lifetime r, be justified in a 
mathematically rigorous way? Only recently has this relation been experimentally tested to 
an accuracy that goes beyond the Weisskopf-Wigner approximation [0,|3|]. 

The other question concerns how the quantum theory of quasistable states can incorporate 
the following notion: there can be no registration of decay products by a detector at times 
prior to the preparation of the decaying state. In other words, we expect that precursor events 
should be assigned zero probability in a theory of quasistable states. We refer to this common- 
sense idea as causality or the preparation- registration arrow of time Mathematical results 
within Hilbert space quantum mechanics challenge the applicability of this notion to quantum 
phenomena []5| and also challenge the validity of the exponential law M. 

In the following two sections, we address the above questions in reverse order and find 
a connection between them, aided by the classic example of Fermi's two atom problem [J7|. 
Then the next two sections present a solution to these questions, defining the Gamow vectors 
for the representation of quasistable states. To make their definition rigorous, and to incor- 
porate simple answers to the above questions, entails selecting separate boundary conditions 
for the space of states and for the space of observables. This requires a slight modification 
of the Hilbert space axiom of standard quantum mechanics. The consequences, physical and 
mathematical, of modifying this axiom so as to incorporate causal decaying state behavior 
are discussed in the conclusion. 

2 Probabilities for Precursor Events 

In his classic review of Dirac's quantum theory of radiation 0j, Fermi postulated a problem to 
test whether the theory satisfied what he considered a sensible requirement, finite- velocity 
signal propagation (or Einstein causality as it is occasionally called). As he posed the 
problem: 
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Let A and B be two atoms; let us suppose at the time t = 0, A is in an excited 
and B in the normal [ground] state. After a certain time A emits its energy which 
may in turn be absorbed by the atom B which then becomes excited. Since the 
light needs a finite time to go from A to B, the excitation of B can take place 
only after the time r/c, r being the distance between the two atoms. @ 

Fermi goes on to solve this problem, making the assumption that the mean life of the excited 
state of A is short and that the mean life of the excited state of B is very long. These facts 
are used to justify several simplifications; in particular he uses them to justify extending the 
range of integration for the emitted photon frequency from [0, oo) to (— oo, oo). Because of 
this approximation, several integrals become exact and he achieves his desired result: there 
is zero probability that atom B will be excited for t < r/c or equivalently for t — r/c < 0. 

As we would say today, Fermi achieved his result by analytically extending the photon 
energy range. Though this extension of energy appears a quantitatively minor modification, 
it has now become clear that precursor events have zero probability only because the integral 
over the physical values of energy < E < oo was extended to an integral over the range 
— oo, < E < oo. The absence of precursor events is an artifact of Fermi's modification of the 
energy range, not a consequence of the quantum mechanics of his time. 

The root of this problem is not that something "propagates faster than light" to give 
non-zero transition probabilities for t < r/c. In fact, if one does not extend the lower bound 
of the energy range to — oo, the transition probability is non-zero even for t < 0, i.e., for times 
before the atom A is excited. This result does not indicate a violation of Einstein causality so 
much as a violation of the basic notion of causality expressed by the preparation-registration 
arrow of time mentioned in the introduction. 

This result is not specific to Fermi's two-atom problem, but is the consequence of a general 
theorem by Hegerfeldt ||: The transition probability between two Hilbert space vectors is 
either identically zero for all times (and therefore, there is no decay at all) or it is different 
from zero for all time — oo < t < oo. This means that excited states cannot be prepared 
at a finite time t > —oo and subsequently decay, although this is the typical situation 
and the case considered by Fermi. This Hilbert space theorem has two consequences: 1. 
non-zero probabilities for precursor events, such as decay products detected at t < t before 
the state was prepared, a result which violates causality, and 2. non-zero probabilities for 
detecting decay products for times t < to + r/c for any finite but arbitrary distance r, or 
equivalently non-zero probabilities for r > c(t — to), which violates Einstein causality. This 
latter consequence was the main concern of || though the former consequence more directly 
conflicts with the concept of causality in general Q, since it does not also involve space 
translation and the constant c of special relativity. Since Einstein causality involves the 
constant c one can only do full justification to this problem in a relativistic theory using 
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Poincare transformations, and not just time translations generated by the Hamiltonian. We 
shall briefly discuss this at the end of this paper and in the Appendix. 

The theorem of || requires that the Hamiltonian be bounded from below, so Fermi 
circumvented the conditions of this theorem by extending the energy values for the spectrum 
of the Hamiltonian H from {E\0 < E < 00} to {E\ — 00 < E < 00} and obtained his result: 
the probability that atom B is excited before t = 0, the time at which excited atom A 
had been originally created, is zero. Within standard quantum mechanics the concept of a 
time to before which the Born probability is zero is precluded because the (Heisenberg or 
Schrodinger) equations of motion integrate under the Hilbert space boundary conditions to 
the unitary time evolution group, valid for —00 < t < 00 ||. 

Fermi's procedure of extending the energy range can be justified as an analytic con- 
tinuation of the energy, but this requires a modification of one of the axioms of quantum 
mechanics, the Hilbert space boundary condition. In fact, Fermi's approximation points the 
way to selecting boundary conditions on the spaces of states and observables in which the 
Schrodinger and Heisenberg equations do not integrate to a group, but to a semigroup. 



3 Cross Section and Exponential 

There are two experimental signatures of quasistable states, the Breit-Wigner or Lorentzian 
for the energy dependence of the decay amplitude 

and the exponential for the time dependence of the decay amplitude 

f(t) = e -iE R t e -r R t/2_ ( 2 ) 

From the point of view of physical observation, the two quantities T in ([!]) and Tr (or 
Tfj/h if the unit of time is not inverse energy) in (|2]) are fundamentally different quantities 
and are measured in different ways. The width T is measured by the Breit-Wigner lineshape 
as, for example, in the cross section of a resonance scattering process: 

The inverse lifetime Tr = 1/t, which is the initial decay rate (considering only one channel) , 
can be measured by fitting the counting rate of decay products to the exponential law: 

i^r K e-Ve- (4) 
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where dN(ti) is the number of decay products registered in the detector during the time 
interval dt around tj. 

The energy scale of the particular quasistable state usually determines which of these 
methods is used to determine the characteristic lifetime or width parameter^ For long-lived 
states, e.g. r > 1CT 8 s (or h/r < 1CT 7 eV), it is reasonably easy to directly measure the 
lifetime by a fit to the exponential and such states are conventionally called decaying states. 
The range of experimentally accessible lifetimes can be extended by exploiting relativistic 
time dilation for fast-moving decaying states; for example, a direct lifetime measurement has 
been made of the decay of the tt° ||1 1|| , for which a value of r = 8.97 ± 0.22 ± 0.17 x 10~ 17 s 
was extracted. This value, which corresponds to a ratio r / 'Er ~ 10~ 7 , is the limit on lifetime 
measurements. For broad resonances, T/Er ph 10 _1 — 10 -4 , the lineshape is easy to measure 
and requires an energy resolution of the detection apparatus (and an energy distribution of 
the prepared quasistable states) comparable to V /Er. They can be resolved for much smaller 
ratios in particular physical systems, such as nuclear resonances, for which widths have been 
measured with the Mossbauer effect for at least T/Er < 10~ 15 (see table 3.1 of |L2] for early 
data). 

Although if one can measure the exponential decay rate, the linewidth will be extremely 
narrow and similarly, although resonances broad enough for their true width to be mea- 
sured have very short lifetimes, there have been quastistable states for which both width 
and lifetime have been measured and the lifetime-width relation T = Tr(= H/t) can be 
experimentally tested. For example, using the Mossbauer effect 40 years ago, the width of 
the first excited state of Fe 57 was found to be 4.7xl0~ 9 eV, which agreed within 10 percent 
with the direct lifetime measurement of r = 1.4 x 10~ 7 s (see section 4.2 of [J12|]). 

A recent precision measurement of V = 9.802(22) MHz = 4.0538(91) x 10 -8 eV for the 
natural linewidth of the 3p 2 P 3 / 2 excited state of Na (using trapped, ultracold atoms f2j) 
made it possible to compare for the same atomic state the lifetime calculated from T, r = 
h/T = 16.237(35) ns, with the lifetime measured directly r = h/Y r = 16.254(22) ns (using 
beam-gas-laser spectroscopy B). The agreement of those two values to an accuracy that 
exceeds the accuracy expected of the Weisskopf-Wigner approximation provides sufficient 
reason to seek a mathematical theory that justifies the relation r = Tr = 1/r as an exact 
equality^. 

1 There are other experimental methods to determine the lifetime of a quasistable state, such as using the 
Primakoff effect for photoproduction of neutral mesons || or mapping the resonant dipole-dipole interaction 
potential for atomic states JlTJ. We will not consider these methods and their theoretical connection to the 
lifetime and width measurements described here. 

2 Note that T = Tr = 1/r has never been verified in the regime of high energy (relativistic) particle 
physics. This is of particular interest for the Z-boson resonance, where there is much debate about the 
parameterization of the lineshape in terms of mass and width and the relation of the width to the inverse 
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Theoretically, the width T and inverse lifetime are often related to each other in 
the following heuristic way: the time evolution of a state vector in standard quantum 
mechanics with a self-adjoint Hamiltonian H is given by 

(j)(t) = e~ iHt ct) for - oo < t < oo. (5) 

To obtain from (JD the exponentials of (§) and @, one takes for the time evolution of a 
decaying state 

(f)(t) = e-^R-iTR^t^ f or _ oo < t < oo (6 a ) 

which would imply that 

H<j> = {E R -iT R /2)<j>. (6b) 

Ignoring the fact that the statements (§) cannot hold for a self-adjoint H in standard quan- 
tum mechanics, one can calculate the probability amplitude as a function of time for finding 
decay products described by if) in the decaying state (f>(t) (|6a| ) as 

0(t)) = e - iHt (j)) = fit) = e- iERt e- rRt/2 (ip } 0), for - oo < t < oo. (7) 

Although it has been known for some time that mathematical theorems of Hilbert space 
quantum mechanics forbid amplitudes with the exact exponential time dependence ||[14 



this heuristic approach is used to justify the exponential counting rate (£D for decaying states. 
However, more is required if the quasistable state is to have both the experimental signatures 
(|) and©. 

To relate the exponential decay amplitude (0) to the Lorentzian, the Fourier transform 
of the vector <p{t) is taken: 

/+oo r+oo 
dt(f)(t)e iuJt = <f) dt e -i{E R -iT R /2)t e iut (g) 

-oo J —oo 

% 



u-{E R - iT R /2) 



Since the vector xi 00 ) is the function (|1|) multiplied by the vector <f>, one argues that the 
relation between ([[]) and (0), and thus the equality T = T R = 1/r, has been established 
theoretically [0]. 

From the point of view of mathematical rigor, there are problems with the above heuris- 
tic derivation: there is the conflict between the self-adjoint Hamiltonian and the complex 

lifetime (see Jl3) and references thereof). 
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eigenvalue in (||) and there is a problem with the boundaries of integration. The Fourier 
transformation is defined by 



+00 

iu>t 



F[f{t)\ = g{u) = / dt /(f)e«* (9a) 
J — 00 

1 /-+00 

T-\g(u)\ = fit) = — dug{u)e-^ (9b) 

^ J-oc 

and this is what is usually used in (§). However, the exact mathematical relation between 
the Lorentzian and exponential is really given |L5| by 



u-{E R - iV R /2) 

and its inverse is 



r+co 

/ dt e -^Rt e -T R t/2 e iut ( 1Qa ) 
JO 



2tt co - (E R - iT R /2) U[t)e \0 for t < ' (1UI) 

The ranges of integration for the variables in the mathematical relation (pT|) are not the 
accepted ranges for the corresponding physical quantities in standard quantum mechanics. 

The energy, the spectrum of the self-adjoint Hamiltonian operator, has a lower bound, 
say Eq = 0, because the scattering energy is a positive quantity and the stability of mat- 
ter requires that the spectrum of the Hamiltonian be bounded from below. However the 
mathematical relation ( |10b| ) requires the energy range —00 < uo < +00. Further, the time 



variable in quantum mechanics extends over — 00 < t < 00 because Hilbert space states 
evolve in time by a one parameter unitary group @ and this time evolution is reversible. In 
the mathematical relation ( |10a| ), the values of t range only from to 00. So the connection 
( [TOP given by the Fourier transform between the Lorentzian energy distribution and the ex- 
ponential time dependence of the wave function requires ranges for time t and energy uj that 
are incongruous with what standard quantum theory allows. 

Tellingly, this energy range —00 < uj < +00 is exactly what Fermi needed for his deriva- 



tion of zero probability for precursor events. Using the relation (|10D, the probability ampli- 
tude is not given by (0), but by 

f (ib d))e~ iERt e~ rRt / 2 for t > 
and the notion of time asymmetry, of an initial time t = before which there is no decay 



probability, emerges from the mathematics (10). 
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In order to define such an extension of the energy range as an analytic continuation of 
the energy, we must require certain analyticity properties of the energy wave functions, and 
that requires modifying the Hilbert space boundary condition hypothesis. This modification 
lead to the definition of the Gamow vector, a mathematical representation of unstable states, 
which we shall use in place of in @. This Gamow vector has both a Breit-Wigner energy 
distribution and an exponential (and asymmetric as indicated by the 9(t) in ( |10b|) ) time 
evolution. 

4 Gamow Vectors 

In this section, we will define the Gamow vectors to be those states whose energy distribu- 
tions and time dependence exactly fulfill the Fourier transform relations (|10|). The decaying 
Gamow vector has a Breit-Wigner energy wave function that extends over — oo < uo < oo 
and is related by Fourier transform to an exponential that starts at finite time to — and 
extends only into the future, t < t < oo, as required by causality. Although these Gamow 
vectors cannot be elements of the Hilbert space, they can be rigorously defined as continuous 
antilinear functionals on a dense subspace of the Hilbert space. They are, like Dirac kets, 
examples of generalized vectors and their use entails the generalization of the Hilbert space 
scalar product. In parallel with our development of the Gamow vector, we will show how 
the same concept of generalized vectors has been used to give meaning to the Dirac kets. 

In standard (time symmetric) quantum theory, a state vector <p is an element of the 
Hilbert space, (f) E7i (and one also assumes a one-to-one correspondence between elements of 
7i and pure physical states, c.f. (^) below). But in Dirac's formalism of quantum mechanics, 
the state vector is expressed in terms of an energy wave function 4>(E) = (E\<f)) by Dirac's 
basis vector expansion 

POO 

4> = / du \u)(u\<p), (12) 
Jo 

where |u;) are the Dirac kets which fulfill 

H\E) = E\E) for E G spectrum(F) = {E\0<E< oo}. (13) 

The spectrum(iZ) is the positive real line 1R + ; it is bounded from below by Eq = 0. 
The scalar product of a vector ip G TC with (f> G TC is 

(V#> = ftM) = / duo (^>H0) = / du^uj)^). (14) 
Jo Jo 
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In other words, the scalar product is given by a (Lebesgue) integral over the physical energy 
values < E < oo in the space L 2 (M + ) of energy wave functions with (E\(j)), (E\ip) e L 2 (IR + ). 
The space L 2 (R + ) is the representation of the physical Hilbert space H. 

The Dirac ket \E) is not actually an element of TC. The Dirac kets can be made mathemat- 
ically rigorous by defining them as continuous antilinear functions (often called functionals) 
on a dense nuclear subspace $ of the Hilbert space, $ C Ji. The Dirac basis vector expansion 
(|12|) is then the Nuclear Spectral Theorem and holds for all <fi £ For example, a possible 
choice for $ (and the most common) is such that <p(E) = (E\<f)) e iS(IR + ), where iS(M + ) 
is the Schwartz space of rapidly decreasing, infinitely differentiable (smooth) functions re- 
stricted to the positive semiaxis 1R + . Then the ket \E) has meaning as an element of $ x , the 
space of continuous antilinear functionals on $, and its delta function energy distribution 
(E\lu) = 5{uj — E) is an element of (<S(1R + )) X , the space of tempered distributions. If and 
ip are elements of $, then their wave functions are smooth and the integral in the scalar 
product ( [Hp is a Riemann integral. If <fi and tp are elements of 7i, the integral in (|TJj) is a 
Lebesgue integral because some elements of the complete Hilbert space L 2 (M + ) will not be 
smooth. 

With the above preparation, we now want to define a vector called ip G (G for Gamow fT6|| ) 
whose energy wave function is a Breit-Wigner distribution, 

iP = /*, WH ^> = fa \u) ^_L_^ . (15) 

We have two alternatives for the boundaries of the integration. If we stick with the rules of 
standard quantum mechanics we have to choose the boundaries of integration according to 
(|12|) to be < uj < oo. If we want the Fourier transform to be the exponential, so that tp G 
has an exponential time evolution then we have to choose the range to cover — oo < uo < oo. 

For values of Yr/Er ~ 1CT 1 and less, there is not much numerical difference between 
these two choices. For the choice < uj < oo, which we call the vector ipappr.i the scalar 
product can be calculated as 

rco poo j-i j 

(^appr.^appr-) = (^appM (^appr.) = I duJ 



2lT (w 



I 




l 






[2E R 








The vector ipap pr * s a vector in 7i, but it is not in the domain of the Hamiltonian H. 
For the choice — oo < uj < oo, which we call ip G , we obtain 

du V 1 ! f ^2 7 (r 7^2 = L ( 16b ) 

^ 2tt (uj - Er) 2 + [Y R /2y 
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Although the numerical difference between ( 16a ) and ( |16b| ) is small, the mathematical dif- 



ference is enormous. Although the quantity (ip G ,ip G ) is a scalar product in L 2 (IR), it is no 
longer a scalar product in the Hilbert space 7i represented by L 2 (M + ), since in ( |Hj) the scalar 
product in that space is defined by integration only over the range < E < oo. 
The difference between the energy distributions is 



r R 1 



K^Wawr.n \ for-oo<£<0 K ' 



and 



\(E\^ u )r = — - -f— — for - oo < E < oo. (17b) 

IV W 71 2tt (E - E R f + (IV2) 2 1 ; 

For small values of Yr/Er the numerical difference between ( |17a| ) and ( |17b| ) for values 
E < is negligible. This is why Fermi may have thought he could neglect this difference 
and make his approximation by extending the lower limit of integration to — oo without 
qualitatively affecting the results. But this can be considered a mistake, because there is an 
important difference between ipa PP r. an d V' G - While the vector ip G ppr ^ * s a vector in the Hilbert 
space ((E\ip G ) e L 2 (M + )), it is not an eigenvector of the Hamiltonian H and it does not 
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have exponential time evolution [ l7| . If one wants an eigenvector of a self-adjoint H with 
exponential time evolution, then the lower limit of the integration must be extended to — oo. 
Therefore, we follow Fermi and extend the boundary of integration in fll5|) to define ifj G : 

?P G = [°du \uo-)(-u)U G ) = Haw \uj-) ( i\— 7 —I . (18a) 

7-oo J-oo ' ' \ V 2tt u - (E R -iT R /2) J 1 ; 

This vector ip G is now no more an element of the physical Hilbert space Ti L 2 (IR + ), but 
ip G is a generalized vector with ideal Breit-Wigner energy distribution that extends over 
f — oo, oo) 




-E\ij G ) = 77-^ — z— - for - oo < E < oo. (18b) 



Unlike the ordinary Dirac kets \lu) in ([T^), we have denoted the kets in (18) by This 
notation has been chosen in conformity with the notation for the solutions of the Lippmann- 
Schwinger equation with — ie in the denominator (the out-going plane wave solutions). We 
will give their exact definition after further developing the notion of a generalized vector. 

The Gamow vector %p G is generalized in the sense that it is not an element of the physical 
Hilbert space 7i, but, like Dirac kets, it has meaning as a continuous antilinear functional 
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(ip\F) = F(if)) on a dense subspace of the Hilbert space, i.e. for -0 G $ C 7i. Generalizing 
the scalar product (i[),4>) = 0(0) to the bra-ket (ip\F) = F{ijj) introduces a larger class of 
vectors F G $ x D 7i x than the continuous antilinear functionals defined by the Hilbert space 
scalar product (-0,0) = 0(VO- This results in a triplet of spaces, a Rigged Hilbert Space, 
$ C H C $ x ; the kinds of generalized vectors \F) in the dual space $ x are determined by 
the choice of $. 

For the Hilbert space vectors 0, the scalar product with all ip G 7i is a functional 
0(0) = (0|0) = ("0;0) on an "0 G 7Y. It is defined in (|14|) (using Lebesgue integration) and 
makes mathematical sense for every ip &7i. However, for the Dirac ket \E) G $ x , 

poo 

\E) = duj \u)S(u-E), (19a) 
Jo 

the functional (V'l-E) on some vector ip makes mathematical sense only if -0 is from a smaller 
space $. For such -0 G $ C 7i, the wave function 0>(-E) 

/"OO 

= = / duj (iP\lu)5(uj - E) (19b) 







is a smooth function in the Schwartz space, (E\ifj) = (ip\E) G iS(M + ). Then, the meaning of 
Dirac's eigenket equation (|T3|) is that of a generalized eigenvector equation 



(Hij)\E) = (ip\H x \E) = E(4)\E) for all ip G $, (20) 

where if x is the unique extension of the operator = H to the set of functionals Fe$ x . 

The vector ip G defined by ( |18a| ) is still more generalized than the Dirac ket \E) and the 
generalization of the scalar product, the bra-ket (ip~\ip G ), makes sense only for every ip~ in 
a still smaller subspace $ + C $ C Ti. The space $ + will be chosen such that the energy 
wave functions (~E\ip~) will be smooth Hardy functions in the upper half complex energy 
plane, i.e. {~E\tp~) G TL 2 + fl <S|r+. Hardy functions are those smooth functions (~£7|^~) 
which can be analytically continued into the upper complex plane, and (ip~\E~) = (~E\ip~) 
can be continued into the lower half complex plane and they vanish at the infinite semicircle 
sufficiently fast |L9||. Only for these 0>~ G does the bra-ket (ip~\ip G ) make sense; i.e. we 



must have ip G G and then the value of the functional ip G at ip is 

w " ' 2iri_ 0o " 1 ' w-(E R -iT R /2) 2m J c o> - (E R - iT R /2) k ' 

In the second integral the closed contour C is from +oo to — oo and then along the infinite 
semicircle; the equality follows because (tp~\uj~) vanishes on the infinite semicircle if ip~ G 
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In distinction to the standard Dirac kets \u) G $ x , the Dirac kets \u>~) used in ( PI] ) (and 
before in ([18])) are functional over the Hardy space $ + , i.e. |u;~) G $*. This we take as the 
mathematical definition of the out-going plane wave solutions of the Lippmann-Schwinger 
equation. From $ + C it follows that <J> X D $ x and so the set of \uj~) is larger than the 
set of standard Dirac kets; the functional (~u\ip~) with i/j G $ + are defined for any complex 
number u of the lower half complex plane. There are also Hardy spaces <3>_ for which the 
kets \uj + ) are defined for any complex number uj of the upper half complex plane; the role 
of these spaces will be discussed in the final section. 



Evaluating formula fl2T|) using the Cauchy formula ||19| , we obtain 

..ox _ * n.. /./.-i. .-x 



m^ G ) = f r&» (v>-> — ,r^ r /9 , = mE R -iv R /2-)^KF R . (22) 

27T J.oo UJ - (E R - li R /2) 



This means that ip G = \/2irT R \E R — iT R /2~) G is something like a Dirac ket extended 
to the complex value E R — iT R /2. Also, comparing this with equation ( |19b| ) shows that the 
Breit-Wigner distribution i/2ir(E — (E R — iT R /2))~ 1 is something like a Dirac 5-function 
when applied to the well-behaved Hardy functions, i.e. just as integrating over 5(u — E) 
maps every function (if)\u) G S to its value at E, integrating over i/2n{uj — (E R — iT R /2))~ l 
maps every function (ip~\u~) G 7i 2 _ PI <S|k+ to its value at E R — iT R /2 with Y R > 0. 

The properties of well-behaved Hardy class vectors ensure that if ip~ G $+, then also 
Hif)~ G $+. Then replacing t/>~ in ( p2|) with if?/' - (or equivalently replacing (ifj~\u)~) with 
(Hip~\uj~) = (tp~\H x \u)~) — u(i/}~\u~)) one proves 

(Htfj-\^ G ) = (^-\H X \^ G ) = (E R - iT R /2)(ip-\il; G ) for all G (23) 

This proves that the Gamow vector ip G is a generalized eigenvector of the self-adjoint and 
semi-bounded Hamiltonian H with a complex eigenvalue (note that this cannot be proved for 
ip G ppr ', it is not an eigenket of H). Omitting the arbitrary ip~ G in (|2"B|), the eigenvector 
property is written in the Dirac notation as 

H*ty G ) =H x \E R -ir R /2-)y^r R ~= {E R -iT R /2)\i) G ). (24) 

Such eigenvectors (eigenkets) of a self-adjoint Hamiltonian with complex eigenvalue do not 
exist in the Hilbert space or even in $ x , the dual to the Schwartz space, but they are in 
D the dual to the well-behaved Hardy class space $+ of the upper half complex 
plane. 
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5 Exponential Decay Law from Breit-Wigner Distri- 
bution 



As shown above, the Gamow vectors are eigenvectors of a self-adjoint Hamiltonian with 
complex eigenvalues, thus justifying the heuristic equation fl6b| ). However this is not enough 
to justify the exponential time evolution (|6a|); for that we return to the Cauchy formula. 

Note that if (ip~\u~) G 7i 2 _ fl <S|r+, then exp(— iujt){ip~ \uj~) will also be in Ti 2 _ D «S|r+, 
though only for t > 0. For t < 0, the exponential factor blows up as uj approaches the 
infinite semicircle. With this fact, the time evolution of the Gamow vector is calculated for 
all ip~ G $+ as follows: 

(e iHt i>-\E R -ir R /2) = ^-\e- mt \E R -iT R /2-) 

"+00 



i r+oo 

— J dw (r\e- iHt \u-) 



oo 
+oo 



LO - (E R - lT R /2) 
1 



= — / duo e- lwt (i)-\uo-) , 

2ttJ_ 00 W 1 >u-{E R -iT R /2) 

= e- iERt e- TRt/2 (^\E R -iT R /2-) (25) 

for all ^ G $ + and for t > only. 

The last step in equality (|25| ) is the generalization of the formula for the Fourier transform of 
the Lorentzian (|10b| ), which one observes if one multiplies the Lorentzian in ( 10rj|) by a Hardy 
class function of the lower complex energy plane, i.e. by a function (ip~\u~) G TC 2 ^ fl <S|r+: 

f /lc^ jf^Q; j- = 9{t)e-< E *-^(r\E R -iT R /2-) (26) 

2vr y_ oc w - {E R - iY R /2) 

for any (?/>~|co>~) G 7i 2 fl <S|r+. 



The mathematical result (ESI) for the Gamow ket ?/> G = \E R — iT R /2 )^/2~nY R can also 
be written in the Dirac notation if one omits the arbitrary tp~ G $+: 

/ + OO • -ILut 

^ 2ir to - [E R - il R /2) 
= e- iERt e- TRt/2 \E R - iT R /2~) for t > only. 

The equality in the second line of (p7|) is the vector analog of the Fourier transformation 
formula ( |10b| ), and fl26|) is the wave function analog of (|10b| ). It expresses the exponential law 
for the time evolution of the vector ip G (t) with idealized Breit-Wigner energy distribution. 
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Thus the Gamow vector is needed to accomplish the task posed: derive the exponential 
law from the Lorentzian energy distribution. The well-known Fourier transformation formula 
( [TOD suggests that we should follow Fermi's "mistake" and extend the energy to — oo. As 
can be expected from the 9(t) in fllPp, time asymmetry about the time t = (which can be 
any time t ^ — oo) results. 

Such a privileged moment of time t = cannot exist in standard quantum mechanics. 
The time evolution in the Hilbert space is necessarily given by the reversible, unitary group 
of operators U^(t) = exp(-iHt) ||; with every U'(t) there exists also an inverse (U^it))^ 1 = 
W(—t). However, in the space the time evolution operator U x (t) = exp(—iH x t) (|27|) is 
the uniquely defined extension of the unitary operator U'(t) for t > only; (U x (t))^ 1 is not 
defined on and U x (t) is a semigroup. The mathematical necessity of this time asymmetry 
and special time t = has a physical consequence: the semigroup time transformation 
properties of the Gamow vector (|27| ) have the notion of time's arrow (or irreversibility) and 
a privileged time to built in. 



6 Conclusion 

In summary, the derivation in introductory textbooks which relates the Breit-Wigner en- 
ergy distribution to the exponential time evolution, and which makes no sense in standard 
mathematical theory of quantum mechanics in the Hilbert space, can be given a mathemat- 
ical meaning. To do so leads to the Gamow vectors and their irreversible semigroup time 
evolution. How one has to construct the Gamow vector is indicated by the mathematical 
relations ( fL0|) for the Fourier transform between Lorentzian and exponential. The relation 
fllOf) already contains the first appearance of time asymmetry and forebodes the inadequacy 
of the unitary time evolution of standard quantum mechanics. The time evolution of the 
Gamow vectors (|27f) can be viewed as just the generalized vector version of the relation ( |10D 
between generalized functions. 

To incorporate generalized vectors requires a partial revision of the mathematical theory 
of quantum physics. The Hilbert space is a vector space with topology (i.e. the notion 
of convergence for infinite sequences of vectors) given by norm convergence. While the 
appropriateness of a linear space of states is experimentally verified by the glorious success 
of the superposition principle, the appropriateness of the norm topology cannot be directly 
tested due to the limited resolution of experimental equipment and the infinite number of 
measurements required to make a topological distinction. It is the algebraic properties of a 
linear scalar product space, and not the topological properties, to which most physicists refer 
when discussing the Hilbert space. The Hilbert space results when this algebraic structure 
is completed with the norm topology, and then the mathematical theory used in |5],|||| is 
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obtained. 

The first revision of this mathematical theory of quantum mechanics was made so that 
elements like the Dirac ket could be incorporated; it is now well-accepted and uses the Rigged 
Hilbert space 

$ C H C $ x (28) 

where $ is usually chosen to be the Schwartz space (which has a stronger topology than TC) 
and then $ x , the space of continuous antilinear functionals on $, is the space of tempered 
distributions . 

The scalar product \(ip,4>)\ 2 in the Hilbert space is interpreted as the probability to 
detect the observable \?p)(ip\ in the state 0. For the generalized "scalar product" with a 
Dirac ket, the physical interpretation of |(i?|'i/')| 2 is as the probability density in energy using 
an apparatus described by ip. Since the macroscopic apparatus will have a smooth energy 
resolution, the amplitude (E\ip) is best described by a smooth function, (E\ip) G <S(R + ) or 
equivalently ip G Thus the mathematic necessity of using the space $ to define the Dirac 
kets has a physical justification. 

The second step of this revision is to incorporate the Gamow vector. In order to include 
generalized vectors with complex (and also negative) energies, the energy wave functions 
[~E\ip~) cannot be just smooth functions but must also be analytically continuable. Instead 
of using the Hilbert space axiom 

{space of prepared states} = {space of detected observables} = 7i (29) 

or the slightly more general revision 

{space of prepared states} = {space of detected observables} = $ C H, (30) 

we distinguish mathematically between states and observables |2l| and make the new hy- 



pothesis [52 



The prepared states are described by: {<fi + } = C C $! (31) 
and the registered observables by: {?A _ } = C 7i C $*. 

Here we use a pair of Rigged Hilbert spaces of Hardy type, where the energy wave functions of 
the vectors <ft + G $_ and ip~ G are well-behaved Hardy functions in the lower and upper 
half complex planes, respectively [^j . The Gamow vectors, together with the out-plane wave 



solutions of the Lippmann-Schwinger equations, are elements of the space $ + . 

The states prepared by the preparation apparatus are, by hypothesis (pT|), not generalized 
vectors like the Gamow vectors (|18D; the states (fi + G G Ti, are represented by the Dirac 
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basis vector expansion (|T^) with wave functions <p + {uj) = that are well-behaved 

Hardy functions in the lower half plane. From (|3l|) one can prove (see the reference in 



citation |19[) that an alternate "complex basis vector" expansion holds for every <p + G 

N r+oo 

+ = El^)^ G ^ + )+ / du\u+)b{u). (32) 
i=i " / -°° 

Here, the ifjf represent N Gamow vectors with eigenvalues z^ = — iTi/2 which are 
associated to N first order resonance poles at the pole positions of the analytically 
continued S- matrix. However, in addition to the sum over the resonance states ipf, there 
appears an integral over the continuous basis vectors with weight function b(u). Whereas 
each Gamow vector in ( |3"2"D corresponds to a Breit-Wigner resonance amplitude, the integral 
corresponds to a slowly varying background in the scattering amplitude. 

Considering just one resonance (N = 1), then the first term of (^) represents the state of 
the resonance with its characteristic exponential time behavior. However, there is always also 
a background term whose time evolution is not exponential and whose energy-dependence 
b(uj) depends on the way in which the state was prepared. This background term is a 
theoretical necessity but in a particular experiment may be unimportant. The influence of 
this background can explain observed deviations from the exponential decay law MM , which 



had been derived mathematically a long time ago || for Hilbert space vectors, and thus 
also for + e C Ti. But the Gamow vectors are elements of and are not in 7i, and 
therefore their time evolution can follow the time-honored exponential law. 

In the spaces and and in the space $ + and (but not in 7i) only semigroup time 
evolution is defined. The semigroup time evolution is the crucial difference between standard 
time symmetric quantum mechanics in the Hilbert space or in a triplet such as (|28|) and a 
quantum theory that includes time asymmetry. By changing the boundary conditions, i.e. the 



space of allowed solutions of the dynamical equations (Schrodinger or Heisenberg) from (J29 
or |3H|) to ([31]), we obtain semigroup evolution. Time asymmetric boundary conditions of time 
symmetric dynamical (differential) equations are nothing new in physics. The cosmological 
arrow of time and the radiation arrow of time are consequences of such a theory. The idea of 



a fundamental arrow of time in the laws of quantum mechanics is also nothing new, c.f. p5 
where it had to be affixed artificially on top of the time symmetric solutions because the 
Hilbert space boundary condition for the quantum mechanical Cauchy problem excluded 
time asymmetry. With the hypothesis fl3"T|), asymmetric time evolution is a consequence of 
the mathematical property of the Hardy space functions |4| (specifically the Paley- Wiener 



theorem for their Fourier transform |T9[) in the same way as reversible time evolution is a 
mathematical consequence of the properties of Hilbert space functions || . 

The Gamow vectors ip G describe resonance states (without background). They are the 
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only mathematical entity which can exactly combine the properties of Lorentzian energy 
distribution and exponential time evolution. The relation between the lifetime r (from the 
counting rate (^)) and the Lorentzian width T (from the cross section ([])) is precisely and 
exactly r = h/T and this lifetime-width relation holds exactly only for the Gamow vectors. 

The Gamow vectors also predict causal probabilities. As a consequence of their semigroup 
time evolution (25.27) the concept of a finite time t >= oo (represented by the semigroup 
time t — 0) is introduced for a quantum system described by ip G , and before that time 
the quantity (e lHt ip~\il} G ) = \e~ tHt ip G ) = \ip G (t)) does not exist. Extending the 
usual probability interpretation of \(ip,<f))\ 2 to the Gamow vector, the generalized scalar 
product |(V'~I' ? /' G (^))| 2 describes the probability to detect the decaying quantum state ip G — 
\Er — iY r/2^)\/2tiY R at the time t with an apparatus described by ip~ . Lets put this into a 
specific physical context for the purpose of example by considering the process where a K° 
is created by ir~p — > K°A and then decays via the channel K° — > 7t + tt~ . The Gamow vector 
ip G represents the decaying state K° and ip~ represents the the decay products 7r + 7r~ which 
the detector registers. Since the decay products cannot be detected before the decaying state 
ip G has been created (or prepared) at an arbitrary but finite time to = 0, the probability 
to detect the decay product \(ip~\tjj G (t))\ 2 = \(e lHt il)~\il) G )f , makes physical sense only 
for t > 0. This means that for the decaying kaon system, we should expect a non-zero 
counting rate for the decay products tt + tt~ only for times after the K° has been created in 
the reaction ir~p — > K°A and leaves the proton target. Therefore the probability for this 
decay (t))\ must be different from zero only for t > to = and this is precisely the 

time asymmetry predicted by (|25|). Thus the semigroup time t = is interpreted as the time 
to at which the creation of the decaying state is completed and the registration of the decay 
products can begin. For the K° system this time is very accurately measurable because the 
K° is created by the strong interaction with a time scale of 10~ 23 s and it decays weakly with 
a time scale of 10 _10 s. 

The result (||) predicts that the probability \(ip~\^ G {t))\ 2 = \(4>~\U x (t)\4) G (t))f to 
detect decay products ip~ from a state ip G which has been prepared at to is zero for t < t . 
Therefore the probability of precursor events is zero, in contrast to the Hilbert space result 
obtained in || which does predict precursor events. This feature is what we mean by 
causality. 

This prediction applies when the detector l^ - )^ - ! is placed right at the position of the 
decaying particle ip° and is in the rest frame of the decaying state. If the decaying state and 
detector are displaced or have relative motion, then time translation alone cannot describe 
the situation and Poincare transformations must be considered. To represent decaying states 



that undergo relativistic transformations, relativistic Gamow vectors were introduced [26 



Unlike the relativistic stable particle states which furnish a unitary (Wigner) representation 
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of the Poincare group relativistic Gamow vectors furnish only a semigroup representa- 
tion of the Poincare transformations into the forward light cone 



x 2 = t 2 - x 2 > 0, t > 0, (33) 

(c.f. Appendix and references thereof). One therefore predicts in this relativistic theory 
space-time translated probability amplitudes 

(V>-|t/ x (l,(t,x))|V G > = (r\e- tiHXt -^M G ) (34) 



only for those space-time translations (t, x) which fulfill 

This is the relativistic analogue of the time asymmetry t > in (p5|). Therefore, decay 
events of the space-time translated state U x (1, (t, x))|-0 G ) are only predicted to occur for 
t > and t 2 > x 2 = r 2 /c 2 or for r/t < c. This means the probability for decay events 
cannot propagate faster than the speed of light and Einstein causality is obeyed by these 
semigroup representations of the Poincare transformations. This confinement to semigroup 



transformations into the forward light cone follows from the Hardy space hypothesis d3l"l). 
Unitary representations of the Poincare group in the Hilbert space are not restricted to the 
forward light cone ( |33"D and therefore do not fulfill Einstein causality. This can be deduced 
from the theorem in reference || which is based on the Hilbert space axiom (^). The 
Hardy space hypothesis predicts probabilities only for semigroup time evolution (for non- 
relativistic) or Poincare semigroup evolution in the forward light cone (for relativistic) and 
as a consequence both the causality conditions "no registration of an observable in a state 
before that state has been prepared" and the Einstein causality condition "no propagation 
of probabilities with a speed faster than light" are fulfilled. 

Fermi's fortuitous (but at the time unjustified) approximation, which was made so that 
causality was maintained, pointed the way to the Gamow vectors, and the mathematical 
relation of the Fourier transform (|T0|) foretold the time asymmetry. 



Appendix: Poincare Transformations of 
Relativistic Resonance States 

To address Einstein causality for the representation of resonances and decaying states, we 
must consider Poincare transformations of a relativistic Gamow vector. We will therefore 
briefly record some definitions and results. Since this subject exceeds the scope of the present 
paper we shall just review the results here and refer to |26| for details. 

The non-relativistic Gamow vector was in some sense a generalization of the Dirac ket to 
complex energies and so we begin by summarizing the representations of relativistic stable 
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states. Stable particles are described by irreducible unitary representation spaces of the 
Poincare group characterized by the invariant mass squared m 2 and by spin j [[27] . The basis 
vectors | [j, m 2 ]pj 3 ) of an irreducible representation space are usually labelled by a component 
of the spin j 3 and the spatial components p of the 4-momentum p = (po, p). Instead of p one 
could equivalently use the spatial components p of the 4- velocity p = p/m = (7, p) = (7, 7V ) 
and then use the 4- velocity eigenkets \[j,m 2 ]pj 3 ) as basis vectors. 

The relativistic Gamow vector \[j, Sr]pj 3 ~) is defined from the resonance pole Sr = (M — 
iT/2) 2 of the S-matrix for a resonance scattering process 1 + 2— > i? 3 + 4 by 

1 f + °° ds 

|[j,s R ]pj 3 -> = — / — — I [7, s]pj 3 ">. (Al) 
2tt J^ OOII s - s R 

Here s = {pi + P2) 2 is the invariant mass squared and the integration extends along the 
lower edge of the real s-axis in the second sheet (the sheet that contains the pole at s#). 
The definition (|A1|) is the relativistic analogue of ( |18a| ). The kets \[j, s]pj 3 ~) are the out- 
plane wave solutions of the Lippmann-Schwinger equation and are functionals over a Hardy 
space $ + (which is slightly different than the space $ + chosen for the non-relativistic case 
but has the same analyticity properties), whereas the usual Wigner basis kets for unitary 
representations are functionals over the Schwartz space. 

The transformation properties of the relativistic Gamow vectors under Poincare trans- 
formations (A, x) (where x = (t,x)) are given by 

U*(A,x)\ij, s R ]pj 3 -> = e-^*- x - v ^^ 3 (H/(A-\p))|[j;s R ]A- 1 p J r) (A2a) 
only for 

t > and t 2 > x 2 . (A2b) 

The matrix D 3 .,^ is the (2j+l)-dimensional representation of the Wigner rotation H^A" 1 , p) = 

L~ 1 (Ap)AL(p) and L(p) is the standard boost, which depends only on p and not the mo- 
mentum p = y/sRp. It is this property that allows construction of the representation [j, sr] 
by analytic continuation of the Lippmann-Schwinger kets to the Gamow kets, 

\\j,s]pj 3 -)^\[j,s R }pj 3 -) (A3) 

in such a way that p remains unaffected and always real. These representations [j, Sr] are 
the "minimally complex" representations. 
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In the limit as the complex invariant mass squared s# becomes real, s# = (M — iT /2) 2 — > 
m 2 , the transformation formula ( |A2a|) looks exactly like the well-known formula for Wigner's 
unitary representations [j, m 2 ] : 

U*(A,x)\\j,s]ph~) = e~ ip ' x ^2 Dj> sjs (W(A~ 1 ,p))\ \j, s] A -1 pjg - ) , (A4) 

where exp(— ip ■ x) = exp(— i r ym{t — v • x)). However, the important difference between ([A3]) 
and Wigner's transformation formula for unitary representations is that Wigner's transfor- 
mation formula holds for the whole Poincare group 

V = {(A, x) : Ae S0(1, 3), x e R 4 } (A5) 

whereas the transformation formulas (|A4|) for the Lippmann-Schwinger kets and ( |A2a| ) for 
the Gamow kets hold only for the orthochronous Lorentz transformations and space-time 
translations into the forward light cone: 

V+ = {(A, x) : A E S0(1, 3), x E R^x 2 = t 2 - x 2 > 0, t > 0}. (A6) 

The formulas ( |A2a| . |A2b|) are the relativistic generalizations of the transformation formula 
(P7|) so as to include all allowed space-time translations, which, as for the time translations 
in (|27|), form only a semigroup V + . 

Specializing ( |A2a| ) to space-time translations (l,x), we obtain from ( |A2a| ) the space-time 
translated probability amplitude in analogy to (|25j): 

(^"(x)|[j,s R ]pj 3 ") = (U{L,x)il)-\\j,SR]pj3~) 

= e-^ M - ir / 2 )^ x - v )^-|[j,s R ]pj 3 -) for all ^ G $+ (A7) 

but only for t > and t 2 > x 2 because the Hardy space $ + does not contain vectors 
U(l,x)\tp~) for which ( |A2b|) is not fulfilled. Introducing the speed of light, the condition for 
the predicted probability ( |A2b|) becomes t > and c 2 t 2 > x 2 or |x|/£ < c. This means the 
probability for decay events is only predicted for the forward light cone and cannot propagate 
faster than the speed of light. Thus Einstein causality is not violated. 
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